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2 ■ ABSTRACT 

> ■ 
r/5 . 

d . Aims. We studied and compared the long-term average hard X-ray (>20keV) spectra of a sample of twelve bright low-mass X-ray binaries 
^ hosting a neutron star (NS). Our sample comprises the six well studied Galactic Z sources and six Atoll sources, four of which are bright 
. ,_( \ ("GX") bulge sources while two are weaker ones in the 2-10 keV range (H 1750-440 and H 1608-55). 

■ Methods. For all the sources of our sample, we analysed available public data and extracted average spectra from the IBIS/ISGRI detector on 
$H ' board INTEGRAL. 

?3 Results. We can describe all the spectral states in terms of the bulk motion Comptonisation scenario. We find evidence that bulk motion is 
always present, its strength is related to the accretion rate and it is suppressed only in the presence of high local luminosity. The two low-dim 
Atoll source spectra are dominated by photons up-scattered presumably due to dynamical and thermal Comptonisation in an optically thin, 
hot plasma. For the first time, we extend the detection of H 1750-440 up to 150 keV. The Z and bright "GX" Atoll source spectra are very 
similar and are dominated by Comptonised blackbody radiation of seed photons, presumably coming from the accretion disc and NS surface, 
in an optically thick cloud with plasma temperature in the range of 2.5-3 keV. Six sources show a hard tail in their average spectrum: Cyg X-2 
(Z), GX 340+0 (Z), GX 17+2 (Z), GX 5-1 (Z), Sco X-l (Z) and GX 13+1 (Atoll). This is the first detection of a hard tail in the X-ray 
spectrum of the peculiar GX 13+1. Using radio data from the literature we find, in all Z sources and bright "GX" Atolls, a systematic positive 
correlation between the X-ray hard tail (40-100 keV) and the radio luminosity. This suggests that hard tails and energetic electrons causing the 
radio emission may have the same origin, most likely the Compton cloud located inside the NS magnetosphere. 
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1. Introduction 

Low-Mass X-ray Binaries (LMXBs) are systems where a com- 
pact object, either a neutron star (NS) or a black hole candi- 
date (BHC), accretes matter via Roche lobe overflow from a 
companion with a mass M < 1M Q . NS LMXBs can be broadly 
classified according to their ti ming and spectral properties 
dHasinger & van der KlislfT989l) . On the basis of this classifi- 
cation, NS LMXBs are divided in Z sources and Atoll sources 
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from the shape of their track in the colour-colour diagram and 
from the different timing behaviour that correlates with the po- 
sition on the tracks. 

The overall spe ctra of Z sources are very soft 
llBarret & Olivel 120021 and references therein) and can be 
described by the sum of a cool (~1 keV) blackbody (BB) and 
a Comptonised emission from an electron plasma ("corona") 
of a few keV. Instead, Atoll sources perform quite dramatic 
spectral changes: when bright, they can have soft spectra 
(similar to Z sourc es) but they switch to low/hard spectra at 
low luminosities. Titarchuk &Shaposhnikov<2005), hereafter 
TS05, implemented a thorough analysis of spectral and 
temporal properties of the Atoll source 4U 1728-34. They 
show that the low/hard spectra at low luminosities can be 
described by the sum of up-scattered spectra related to the 
Comptonisation of the disc and NS surface soft photons. They 
found that the Compton cloud electron temperature is of the 
order of a few tens of keV. These spectra are very similar to the 
hard-state spectra of BHCs but they are softer (BHC photon 
index r ~ 1.6 ± 0.1 vs NS T ~ 2.1 ± 0.1) as expecte d from 
the theory (see, for example, fntarchuk & Fioritol2004 . TS05 
found that the high luminosity state spectrum of 4U 1728-34 
consists of the sum of two pure blackbody-like spectra with 
colour temperatures of about 1 keV and 2.2 keV. The softer 
BB component is presumably related to disc emission as the 
harder one is related to the NS emission. It is worth noting 
that TS05 also found that when the source undergoes hard-soft 
transition, all power spectrum (PDS) frequencies (QPO and 
break frequencies) increase with the photon index, with no sign 
of saturation. Note that the index-QPO correlation observed 
in BHC show s the index saturation at high v alues of QPO 
frequency (see Sh aposhnikov & Titarchukl20 06). 

So far, only Atoll sources (and more generally X-ray 
bursters) have been observed with low/hard spectra (i.e. 
Comptonising corona of few tens of keV). Z sources always 
have soft Comptonisation spectra (Comptonising corona of 
few keV) and can have an additional hard X-ray compo- 
nent dominating the spectrum above ^30keV. This compo- 
nent is on top of the soft spectrum and is highly variable 
with most of the emiss ion remaining soft ("see barrel 1200 it 
iDi Salvo & Stelfall2002l for a review on NS LMXB spectra). 
Hence, we would expect that the average high energy spec- 
tra of Atoll sources have a strong component above 30keV 
and that in the soft Z sources this component is less promi- 
nent and smeared out in the time averaged spectrum. Z sources 
spend most of their time in the high/soft state, but they may 
show the transition to harder states at lower luminosity. The 
sensitivity of the past missions may have introduced an obser- 
vational bias, similarly to the lack of a continuous coverage 
of the Galactic plane and Centre in the less explored hard X- 
ray range (above 20keV). Moreover, the concentration of these 
sources towards the Galactic centre makes it difficult to observe 
them with non-imaging instruments. Consequently, data anal- 
ysis and interpretation of such observations is extremely prob- 
lematic. All these instrumental biases can be minimised with 
the use of the recently launched INTErnationa l Gamma-Ray 
Astrophysics Laboratory, INTE GRAL dWinkler et alJ l2003l) . 
The imager INTEGRAL/IBIS JUbertini et alJ 120031) has high 



Table 1. LMXBs s tudied in this pape r . D: d istance (in kpc) 
from references in Midiari & Fender (2006) exce pt for (*) 
from I Christian & Swaajd (m& and (**) from | Fordetal. | 
(2000); Rate: average 22^40 keV counts/sec of the source as 
obtained from the mosaic image shown in Fig. ^ Multiply 
by ~10 to obtain a flux estimate in units of mCrab. 
F(lmCrab) 22 -4ofc e y ^.SxlO^^ergs" 1 cm" 2 ; SNR: signal 
to noise ratio in the 22-40 keV band; MaxEn: maximum energy 
channel (keV) with a signal to noise ratio higher than three in 
the average spectrum; T exp : effective exposure time in ksec. 



Source 


D 


Rate 


SNR 


MaxEn 






(kpc) 


(cps) 




(keV) 


(ksec) 


Z sources 


Sco X-l 


2.8 


58.3 


1757 


150 


266 


GX 340+0 


11 


2.4 


124 


46 


433 


GX 349+2 


5 


3.4 


189 


39 


265 


GX5-1 


9.2 


3.8 


330 


80 


1091 


GX 17+2 


14 


4.4 


203 


80 


248 


Cyg X-2 


13.3 


2.2 


57 


55 


149 


Atoll sources 


H 1608-522 


4** 


0.8 


36 


150 


299 


H 1705^140 


i j** 


3 


144 


150 


307 


GX9+9 


5* 


0.97 


55 


35 


139 


GX3+1 


5.6* 


0.96 


86 


43 


2027 


GX9+1 


7* 


1.2 


91 


37 


491 


GX 13+1 


7 


0.96 


57 


80 


290 



sensitivity, about ~10 times better than G/M/VA7YSIGMA, 
coupled to imaging capability with 12' angular resolution 
above 20 ke V. 

In this paper we report the study of the average hard X-ray 
spectra of twelve NS LMXBs perfor med with the low energy 
(20-200 keV) IBIS detector, ISGRI jLebrun et al.l2003h . using 
a coherent and large sample of data, free from systematic ef- 
fects which play a role when combining data from different 
missions. The sample of the LMXBs chosen is given in Tabled 
and comprises six Galactic Z sources and six Atoll sources, 
four of which are bright ("GX") bulge sources while two are 
weaker ones in the 2-10keV range. 

Our approach is two-fold: on one side, for comparison pur- 
poses, we study the average spectra in terms of phenomenolog- 
ical models as done in the literature, on the other, we study the 
sources in the frame of a physical model in the attempt to find 
a self-consistent scenario that describes all the spectral proper- 
ties we observe. We discuss the similarities of such a scenario 
with the black hole LMXB case as well as the radio - X-ray 
correlation that is typical of LMXBs. 

2. Observations and data analysis 

We have analysed INTEGRAL data publically available in 
which the sources in Table [2 were in the Fully Coded Field 
of View (FCFOV) of IBIS. The additional criterion of a mini- 
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ISGRI 22-40 keV 




(-5 Msec) 



* ♦ 



10.000* 0.000 350.000 340.000 33 

Fig. 1. IBIS/ISGRI 22— 40keV mosaic of the Galactic bulge (about 5 Msec). Only the sources studied in this paper are labeled. 
The location of Cyg X-2 and Sco X-l is not covered by this map. 



mum of good time of 1000 sec for IBIS/ISGRI led to a total of 
2263 pointings each with variable exposure time (from about 
1800 sec up to about 3600 sec) spanning from January 2003 
to May 2004. Version 5.1 of the Off-line Scientific Analysis 
(OSA) software has been used to analyse the data. The descrip- 
tion of the algorith ms used in the IBIS/ISGRI scientific analy- 
sis can be found in lGoldwurm et all J2003I) . 
For each pointing we extracted images in the 22-40 and 
40-80keV energy bands. The images were used to build light- 
curves as well as a final mosaic. The mosaic images revealed 
the sources that were active in the field of view at the time of the 
observations. The list of all the detected sources (one per source 
of interest) was then used in the spectral step, where spectra for 
all the active sources in the field of view were simultaneously 
extracted with the standard spectral extraction method. The sin- 
gle pointing spectra were then averaged into one final spectrum 
per source. These final average spectra have been used in our 
spectral study adding 1 % systematic error (with the exception 
of the bright Sco X-l for which we added 1.5% systematic 
error). In the fit, only data points with more than 3<r in the 
22-200 keV range have been considered. 
To cross-check our results, we also used the alternative method 
of extracting spectra from the mosaics. To do so, we re-ran the 
imaging step in twelve different bins (instead of the 22-40 and 
40-80 keV previously mentioned) for all the pointings and then 
extracted a spectrum using the flux from each energy map. We 
verified that the two different methods give compatible results. 
In this work we show only the results of the former method (the 
standard extraction) for which we have a finer binning. 

3. Results 

Figure ^ shows the final sky map of the pointings around the 
Galactic bulge in the 22-40 keV band. The sources studied 
in this paper are labeled, with the exception of Cyg X-2 and 
Sco X-l that lie outside the image. The 22^-0 keV average 
count rate and detection significance of each source, taken from 



the mosaic image shown in Fig.^ are reported in Tabled The 
count rate is in counts/sec and multiplied by ~10 gives an esti- 
mate in mCrab 1 . 

Figures and show the IBIS/ISGRI (22^10 keV and 
40-80 keV) light-curves as obtained from the imaging analy- 
sis (re-binned to ten day data bins). The simultaneous quick- 
look results provided by the RXTE/ASM team are also shown 
(one day data bins). One Crab is approximately 75 counts/sec 
in the RXTE/ASM 2-12 keV energy band, 93 counts/sec and 
61 counts/sec in the IBIS/ISGRI 22^10 and 40-80 keV band, 
respectively. 

3.1. A phenomenological approach 

As a first attempt we tried to fit all the spectra with the 
same simple model, in order to compare them directly. It 
soon turned out that this was not possible since the long- 
term average spectra identified three main spectral states. 
Phenomenologically, we can classify them in terms of a 
wid ely used the rmal Comptonisation model, COMPTT by 
iTitarchukl J 1994 . We found that the spectra of GX 349+2, 
GX 9+1, GX 9+9, and GX 3+1 are well described by a 
single COMPTT component while we need two components, 
COMPTT+PL (thermal Comptonisation plus power-law), to de- 
scribe the data for Sco X-l, GX 5-1, GX 340+0, GX 17+2, 
Cyg X-2 and GX 13+1. On the other hand, the spectral shapes 
of H 1608-522 and H 1705-440 can be well fitted by a simple 
power-law, PL. 

By comparison with black hole LMXBs, hereafter we call 
the spectral states fitted with a single COMPTT a very soft 
state for which the spectral shape is well described by the 
blackbody-like shape slightly modified by Comptonisation. 
Using this analogy, we call our COMPTT+PL-state an inter- 
mediate state. In fact, the intermediate state in the BHCs is 
characterized by the blackbody-like component at low ener- 

1 The average Crab count rate in different positions of the FCFOV 
in the same band is about 93 counts/sec. 
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Fig. 2. RXTEI 'ASM and IBIS/ISGRI light-curves (counts/sec) of the Z sources. One Crab is approximately 75 counts/sec in 
RXTEI ASM (2-12keV), 93 counts/sec in IBIS/ISGRI 22^10 keV and 61 counts/sec in IBIS/ISGRI 40-80 keV. 



gies and the steep power-law component at higher energies. 
In the low/hard state spectra of BHCs, the blackbody signa- 
ture is smeared out but the power-law is prominent. Keeping 
in mind this observational fact for the BHCs, we can call our 
PL state a low/hard state. The three states are shown in Fig.|4] 
In the plot we hav e also included the spectrum of the Atoll 
source GX 354-0 dFalangaetalJl2006l) to point out that the 
low/hard state of our classification, associated to H 1608-522 
in the plot, is not a "pure" low/hard state as in the case 
of GX 354-0 (thermal Comptonisation with cut-off) but is a 
low/hard state where an additional physical mechanism pro- 



ducing a non-attenuated power-law starts to be important (see 
Sect. 4.1). Nevertheless, for simplicity hereafter we will refer 
to the spectra of H 1608-522 and H 1705-440 as to low/hard 
states. Table 13 summarises the spectral models used to fit the 
data while the corresponding plots with the spectra and best fit 
models are shown in Figs.|5]and|6] 

3.1 .1 . The very soft spectral state 

Three bright Atoll sources (GX 3+1, GX 9+1 and GX 9+9) 
and the Z source GX 349+2 can be well described by a sim- 
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Fig. 3. RXTEI ASM and IBIS/ISGRI light-curves (counts/sec) of the Atoll sources. One Crab is approximately 75 counts/sec in 
RXTEI ASM (2-12keV), 93 counts/sec in IBIS/ISGRI 22-40 keV and 61 counts/sec in IBIS/ISGRI 40-80keV. 



pie COMPTT component. A hard tail has been observed to 
dominate the spec tra above 30keV in a pointed observation 
of GX 349+2 JDi Salvo etalJkOO lb . but it is likely not fre- 
quent/bright enough to be detected in a 265 ksec average spec- 
trum. Indeed, no hard tai l was detected in the spectrum of 
GX 349+2 bv llaria et ail J2004) using BeppoSAX data even 
though the source was in the same position in the Z-track where 
the tail was first discovered. No hard tails have ever been de- 
tected in the spectra of GX 9+9, GX 9+1 and GX 3+1. 

Only recent results have app eared as far as the X -ray spec- 
trum of GX 9+9 is concerned. Kon g et al based on 
RXTE data, report a spectrum that can be described by the 



typical two component model used for bright LMXBs host- 
ing a neutron star: a softer component (e.g. blackbody) plus 
a Comptonised harder one (e.g. COMPTT). We are not sensitive 
to the soft blackbody component with IBIS/ISGRI and our re- 
sults as far as the COMP TT part are concern ed are in agreement 
with what was found bv lKong et ai1 l2006) i.e. a Comptonising 
plasma kT e of about 4keV with a high optical dept h r about 7. 
Our R XTE-INTEGRAL simultaneous observations dVilhu et alJ 
2006| show that the broad band spectrum of GX 9+9 can be 
also well described by the sum of a soft component (DISKBB in 
XSPEC terminology) and thermal Comptonisation (COMPBB) 
by a hot optically thin electron plasma. 
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Table 2. Summary of the spectral properties of the sources. No error means the parameter was fixed to the indicated value. Model: 
best fit model. See text for details; T: photon index of the power-law; kT see( i: temperature of seed photons for Comptonisation in 
keV; kT e : Comptonising plasma temperature in keV; r: optical depth of the Comptonising plasma; Flux^o-io keV) '■ A ux obtained 
between 2CM-0 keV. For the sources detected up to 150keV also the 40-150keV flux is given. 



Source 


Model 


r 


kT seed 

(keV) 


kT e 
(keV) 


T 


X 2 /dof 


F1UX(20— 40 keV) 

(ergs -1 cm -2 ) 


FlllX^4o_l50 keV) 

( erg s -1 cm -2 ) 


Sco X— 1 
GX 340+0 


rOMPTT-t-PT 

COMPTT+PL 


o 85 +o.oi 
[2.5] 


N 31 
[0.93] 


2 87 +0 - 01 

Z - ' -0.02 

[3.0] 


9 45+O-H 

M - 40 -0.34 

7.73+2-?! 

^— 0.78 


42/34 
24/26 


5.8 x 10 -9 
2.4xl0 -10 


2.2x 10 -10 


GX 349+2 


COMPTT 




[1.3] 


9 7 c+0.33 
/ -'°-0.15 


12+ 16 


22/24 


3.5xl0 -10 




GX5-1 
GX 17+2 


COMPTT+PL 
COMPTT+PL 


4 7 +0.6 
' -0.11 

[2.7] 


[1.0] 
[0.6] 


2 42 +0 - 41 

Z '^ Z -0.06 

2.94+0.04 


[11] 


35/29 
35/31 


3.5xl0 -10 
4.5 xlO -10 




CygX-2 


COMPTT+PL 


[1.96] 


[1.09] 


o nn +0.10 
aM£ -0.10 


10.68 


29/27 


2.1xl0 -10 




H 1608-522 


PL 


2.62 ±0.11 








45/36 


8.7xl0 -11 


9.1xl0 -11 


H 1705-440 


COMPPS 
PL 


3.15 + 0.04 


9 70+1-14 
-0.25 


122 tH 


0.18j£S 


40/34 
36/37 


3xl0 -10 


1.9xl0 -10 


GX9+9 


COMPPS 
COMPTT 




2 56 +0 ' 62 

Z ' 0D -0.47 

[1.3] 


47 ±7 
3.66 + 2 5 


0.63:°;" 
6(>0.4) 


31/35 
27/23 


9.2xl0 -11 




GX3+1 
GX9+1 


COMPTT 
COMPTT 




[1.3] 
[1.3] 


o as +0.59 
Z - 88 -0.28 

2 30 +0 - 45 


7 5+ 4 - 7 

'■° -2.9 

26 (> 7) 


46/24 
32/22 


9.3xl0 -11 
1.3xl0 -10 




GX 13+1 


COMPTT+PL 


[2.8] 


[1.0] 


9 7 r+0.25 
Z -'°-0.23 


[8] 


23/22 


8.2xl0 -11 





Table 3. Best fit parameters for the bulk motion Comptonisation model (BMC). No error means the parameter was fixed to the 
indicated value. kTbb'- BB colour temperature, in keV; a: energy spectral index {T=a +1); log A: covering of the BB by the 
Compton cloud ; Flux^o-40 keV)'- A ux obtained between 2CM10 keV. For the sources detected up to 150 keV also the 40-150keV 
flux is given. 



Source 


(keV) 


a 


log A 


X 2 /dof 


FluX(20-40 keV) 

(ergs -1 cm -2 ) 


F1UX( 40 _150 keV) 

(ergs -1 cm -2 ) 


Sco X-l 


z - 3o -0.01 


t 99+0.06 
J - z '-0.02 


-1 47+0.01 
-i.+ '_o.oi 


68/35 


5.8xl0 -9 


2.2xlO -10 


GX 340+0 


2 43+ - 18 

Z - HO -0.31 


<6 


1 19+O.8I 
_1 - lz -1.05 


24/25 


2.45 xlO -10 




GX 349+2 


2 47 +0.04 
z -^'-0.12 


<6 


-1 62+ ' 69 

1 - oz -1.03 


20/23 


3.5xlO -10 




GX 5-1 


7 -ic+0.07 
^.JJ-0.09 


3 84 +0 ' 72 


-1 06+ - 28 


34/30 


3.75xl0 -10 




GX 17+2 


Z - 05 -0.06 


20 +1.46 
z - u -1.41 


1 O9+0.55 


32/30 


4.5xl0 -10 




Cyg X-2 


9 7Q+0.08 
Z - /8 -0.07 


<1.8 


-0.561^ 


26/26 


2.1xl0 -10 




H 1608-522 


l.49±lf 3 


1 50+ - 14 

1- J "-0.16 


-1.24+ 1 ' 45 


39/34 


8.5xl0 -11 


9.6xl0 -11 


H 1705^140 


0.84i°il 


9 , <r+0.05 
z - iJ -0.05 


9 ,, + 10.13 
-Z.1J_ 3 gs 


36/35 


2.9xl0 -10 


1.9xl0 -10 


GX9+9 


7 CT+0.18 
Z - J: '-0.17 


<6 


-0.70t°di 


24/22 


9.4xl0 -11 




GX3+1 


2 30+ 11 

Z - JU -0.14 


<6 


, 97+0.37 
_1 - z '-1.32 


42/23 


9.3xl0 -11 




GX9+1 


9 9n +0.05 
z - zu -0.07 


<3.2 


-1 09+ 016 

luy -0.33 


24/21 


1.32X10 -10 




GX 13+1 


9 ,9+0.27 
z - 3z -0.36 


2.02 + l; 48 3 


1 o 4 +0.58 


23/23 


8.2xl0 -11 





The fit we obtain for GX 3+1 is not good (x 2 /dof=46/24); 
however, as can be seen in Fig. |S] this is not due to a clear 
excess or systematics, more likely it is due to source spectral 
variability. Our average spectrum parameters show that also 
GX 3+1 has a Z - like s pectrum, similarly to what found by 
lOosterbroek et alJ d20Q lh using BeppoSAX data (optical depth 
r about 6 and Comptonising plasma temperature kT e about 
2.7 keV). 



Ilaria et all J2005I) presented the 0.12-18keV spectrum of 
GX 9+1 from a long BeppoSAX observation. The 0.12-18keV 
spectrum could be well fitted by a blackbody model plus 
Comptonised component, the latter with an electron temper- 
ature kT e of 2.3 keV and optical depth t of about 15. No firm 
conclusion could be drawn above 20keV due to source con- 
fusion issues. The INTEGRAL observation above 20keV does 
not suffer from source confusion and shows that the spectrum 
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Fig. 4. IBIS/1SGRI average spectra of GX 354-0 ("pure" 
low/hard state, see text), H 1608-522 (low/hard state), GX 5-1 
(intermediate state) and GX 3+1 (very soft state). 

of GX 9+1 is very similar to the bright Atoll sources GX 9+9 
and GX 3+1, i.e. very similar to the typical Z source spectrum. 
There is a hint of excess in the spectrum of GX 9+1, Fig. [6] but 
it is not significant with the current exposure time. 

3.1 .2. The intermediate spectral state 

A simple COMPTT model was not good enough to describe 
the data due to a systematic excess above ^30keV for 
ScoX-1 (x 2 /dof>60), CygX-2 (45/27), GX 17+2 (61/31), 
GX 5-1 (96/31), and GX 13+1 (38/24). With the addi- 
tion of a power-law, the fits improved significantly (see 
Table |2j with the temperature of the Comptonising plasma, 
kT e , stabilising itself around 2.5-3 keV. Within this sam- 
ple, transient hard t ails had already been observed in the 
spectra of Sco X-l dD' Amico et all l200ll) based on RXTE 
data, GX 5 -1 dPaizis et alJl2005l) based on IN TEGRAL data 
GX 1 7+2 jFarinelli et al.ll2005h and Cyg X-2 JPi Salvo et ail 
l2002h based on BeppoSAX data. In our study, the seed pho- 
ton temperature could not be constrained and was fixed to the 
value [1.0-1.3] keV for all the sources, with t he exception of 
GX 3 40+0 for which it was fixed at 0. 93 keV dLavagetto etail 
l2004 . and GX 17+2 fixed at 0.6 keV dFarinelli et all2005h . 

The spectrum of GX 340+0 could be fitted with a sim- 
ple COMPTT model (y 2 / dof = 32/29) with k T SPPd fixed 
at 0.93 keV as found by Lavag etto et ~M |2004), providing 
kT e =6. 5+2:1 keV and T ~ 1 -' 7 -i.l- We note however, that if we 
freeze kT e to 3 k eV, as found in the Be ppoSAX broad-band fit 
(0.1-100keV) bv lLavagetto et all i2004 . around -40keV the 
residuals systematically deviate from COMPTT and the fit is no 
longer acceptable (x 2 /dof=58/27). In this case the addition of 
a power-law (T frozen at 2.5, as found in BeppoSAX data) is 
required. 

The Comptonising plasma parameters obtained in the fit of 
the spectrum of Sco X-l show values typical of Z sources (kT e 
of a few keV and optical depth between 5-15). The excess left 
by this model is well fitted by a power-law with no measured 
cut-off up to 150keV. Such an extended power-law can be an 
indication for bulk motion (dynamical) Comptonisation or non- 



thermal electron population given the bad fit obtained with 
differe nt thermal Comptonisation models. See iDi Salvo et alJ 
(2006) for a study of the spectral variability of Sco X-l with 
IN TEGRAL. 

Farin elli et alJ d2005l) used the same model used here 
(COMPTT plus power-law) to describe the spectrum of 
GX 17+2, so a more quantitative comparison is possible. We 
fixed the power-law slope (T — 2.7), optical depth (t = 11) 
and seed phot on temperature W T see( j — 0.6 keV) to what was 
found by Farinelli et al. (2005) and find a similar electron gas 
temperature but a lower power-law normalization at 1 keV (our 
~0.23 versus their ~2.1). This is compatible with a picture of 
a fairly constant soft component and transient hard tail compo- 
nent that is at maxi mum when observed in a pointed observa- 
tion (Fari nelli et al.l2005l) and at a lower value when averaged 
on long-term observations (due to its transient nature). 

Similarly, for Cyg X-2 the power-law slope could not 
be constrained and we used the average value found by 
IPi Salvo et alJ d2002h who also used the COMPTT plus power- 
law model. The hard tail was the strongest in the horizontal 
branch and we used the photon index average value from the 
upper and lower horizontal branch (~2). The same was done 
for the seed photon temperature (~1 keV) and optical depth 
(~11). Similarly to GX 17+2, w e obtain an electron ga s tem- 
perature close to what found by Pi S alvo et alJ (12002h and a 
power-law normalization of ^0.02 (to be compared to their 
~0.03). The difference in the power-law normalization is much 
less than for the case of GX 17+2. This could be due to the fact 
that Cyg X-2 spends more time than GX 17+2 on the horizon- 
tal branch and so the average spectrum of Cyg X-2 is closer to 
the horizontal branch spectrum than in the case of GX 17+2. 

The hard X-ray spectrum of GX 5-1 has always been diffi- 
cult to study due to the presence of the nearby bla ck hole candi- 
date GRS 1758-258 (40' apartl. lAsai etaP d 19941) first detected 
a possible hard tail in the spectrum of GX 5-1 although con- 
tamination from the nearby black hole could not be excluded. 
The first non-contaminated observation of a systematic harden- 
ing of the source spec trum in the horizontal branch was done 
bv IPaizis et al J d2005h with INTEGRAL data and the presence 
of the hard tail is confirmed here, based on a larger data set. 

The addition of the power-law component in the spec- 
trum of GX 13+1 required some parameters to be frozen. 
We fixed the seed photon temperature and optical depth to 
the values found by Ho man et al 1 d2004h (1 and 8keV, respec- 
tively). The photon index was found to be steady around 2.8 
(even with a cut-off PL plus PL model and bremsstrahlung plus 
PL). Fixing the slope to 2.8, we obtained the Comptonising 
plasma temperature of about 2.75 keV given in Table [2] The 
natu re of GX 13+1 has alwa ys been a matter of debate: 
iHasinger & van der Klisl Jl989) classified it as a bright Atoll 
source (hence our classification), although they noted that it 
showed properties which w ere closest to those seen in the Z 
sources; Schu lz et alJ (11989) grouped it with the high luminos- 
ity sources that included th e six s ource s that were later clas- 
sified as Z sources; Homan et al. d!998l) discovered quasi pe- 
riodic oscillations (QPOs) similar to the horizontal branch os- 
cillations found in Z sources; based on simultaneous radio and 
X-ray observations. iHoman et al.l ( 120041) showed that the X-ray 
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and radio luminosities of GX 13+1 are actually in the range ex- 
pected for Z sources unlike the Atoll sources GX 9+1, GX 9+9 
and GX 3+1 that have similar X-ray luminosities but much 
lower radio ones. Our work shows that indeed GX 13+1 be- 
haves like a Z source not only due to its variability and ra- 
dio properties but also from the spectral point of view with 
an average spectrum very similar to Z sources, i.e. thermal 
Comptonisation plus hard tail. 

3.1 .3. The low/hard spectral state 

The two remaining Atoll sources, H 1608-522 and H 1705-44, 
are the only ones of our sample that can be well described 
by a single power-law. The burster H 1608-5 2 has been de- 
tected beyond lOOkeV also in the past JZhang et alJ 1 19961 
with BATSE observations) when it was in a very bright state 
(about lOOmCrab, 20-100 keV) with a spectral break or cut- 
off around 65 keV. 

H 1705-44 has been observed t o show dramatic spec - 
tral changes in its X-ray spectrum ( Barret & Olive 2002). 
The 2-80 keV could be well fitted by a blackbody plus 
Comptonisation model (COMPTT) with a temperature of about 
15keV i.e. the spectrum was not detected above 100 keV. In 
this study we detect the source up to about 150keV, that, to 
the best of our knowledge, is a new result. The high sen- 
sitivity of IBIS/ISGRI combined with the long observation 
time (307 ksecs) are essential ingredients for this new detec- 
tion, similarly to wha t happened e.g. in t he case of Aql X- l 
llRodriguez et all2006l) and 4U 1850-087 dSidoli et all2 006). 

To have a feeling of the plasma temperature in 
H 1608-522 and H 1705-440, we fitted their spectra with 
the Comptonisation model by IPoutanen & Svenssonl (|1996) 
(COMPPS in XSPEC terminology) obtaining Comptonisation 
of soft seed photons (~2keV) by an optically thin (t<1) hot 
plasma. We assumed a spherical geometry and a blackbody 
seed photon spectrum (Table|2ji. 

3.2. Generic Comptonisation Model 

As we have mentioned above, in our analysis of X-ray spec- 
tra from NS LMXBs we found similarities with black hole 
LMXBs. The BHC and NS spectral shapes are generic: they 
consist of blackbody-like (BB) and power-law components. 
Sometimes one needs an exponential rollover in order to termi- 
nate the power-law component at high energies. The main dif- 
ference is that NS spectra are usually softer for the same state. 
Another difference in these spectra is that in the NS case there 
are two BB components (not one as in BHC spectra) which 
can be related to the emission from the disc (BB colour tem- 
perature about 1 keV) and NS surface (BB colour temperature 

ab out 2.5 keV), (TS 05). 

iTitarchuketaD i 19961 Il997l) . hereafter TMK96 and 
TMK97, introduced the Generic Comptonisation model 
(BMC model in XSPEC) that takes into account the dynami- 
cal Comptonisation (converging inflow, expected in the vicinity 
of the central object) along with thermal Comptonisation. It is 
worth noting that this model can be applied to fit the observed 



spectra of NS and BHC sources. The BMC model reads as: 

F(E) = ^_(BB + A x BB *G) (1) 

where G(E, Eq) is the Green's function obtained as a solution 
of radiative diffusion equation in which both effects of thermal 
and dynamical Comptonisation are included. BB stands for the 
injected blackbody-like spectrum, namely BB oc B u (T co i), 
BB * G is a convolution of BB with the (Comptonisation) 
Green's function G(E, Eq), Cn is the BB normalization coef- 
ficient (that depends on the compact object mass and distance 
to the source). The factor 1/(1 + A) is the fraction of the seed 
photon radiation directly seen by the Earth observer, whereas 
the factor A/(l + A) is the fraction of the seed photon radiation 
up-scattered by the Compton cloud. 

G(E, Eo) is the response of the Compton medium (cloud) 
to the injection of monochromatic line with energy Eq - TMK97 
showed that G(E, Eq) is a broken power-law with blue wing 
index a and red wing index always more than a + 3. Note that 
G(E, Eq) approaches Delta-function 5(E — Eq) for a»l. 
High values of a are an indication of low efficiency of 
Comptonisation that usually occurs when the plasma of the 
Compton cloud is very close to equilibrium with the seed pho- 
ton environment. The electron temperature of the Compton 
cloud differs by a factor of a few from the BB seed photon tem- 
perature. We also have to note that a single BB spectrum can 
be reproduced either with A <C 1 (small coverage of the seed 
photon radiation by the Compton cloud) or with A ^> 1 (large 
coverage) and a » 1 (inefficient Comptonisation). The free 
parameters of the BMC model (apart from the normalization 
Cn) are the BB colour temperature, kT^, the spectral index a 
(photon index T = a + 1) and log A. 

We can describe all the spectra studied in this paper in 
terms of bulk motion Comptonisation. The results of the fits 
of our INTEGRAL spectra with the BMC model are reported 
in Tabled For GX 340+0, GX 349+2, GX 9+9 and GX 3+1 
only kTbi, is well constrained. This is not surprising given that 
in their spectra only the thermal bump is unambiguously de- 
tected. The best-fit colour temperatures are about 2.5 keV using 
COMPTT or BB models. GX 17+2 and GX 13+1 have simi- 
lar best-fit T-values, even if not well constrained, whereas the 
slope of GX 5-1 is steeper but much better constrained than 
in GX 17+2 and GX 13+1. The upper value of T reported 
for Cyg X-2 (< 3) could be indicative of the presence of pho- 
ton up-scattering (Comptonisation) of BB photons originated 
in the disc and NS surface. One can arrive to similar conclu- 
sions for GX 9+1 even if the effect seems to be less prominent 
there (T < 4). For the two dim Atoll sources H 1608-522 and 
H 1705^440 we find that T is well constrained while there is 
no real information on kT^. This is not surprising given that 
in their spectra we do not see evidence of the thermal bump 
(described by the BB in the BMC model) and so it is impos- 
sible to infer either the BB colour-temperature or its relative 
contribution to the merging flux. In Sco X-l, the whole set of 
the best-fit parameters is well constrained, given that both ther- 
mal bump and the hard tail are detected with high statistical 
significance. 
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Fig. 5. Z sources: average IBIS/ISGRI spectra and best fit models shown in Table|2] 



Since our analysis uses data starting from 20keV, careful 
interpretation should be given to the origin of the BB seed pho- 
tons found in Table H^kT^^ 2.5 keV). A broad band analysis 
of GX 17+2 and Cyg X-2 spectra (Farinelli, private commu- 
nication) seems to suggest that, similarly to the BHC case, the 
seed photons for the BMC model are rather soft, w ith a temper- 
ature close to disk emission values, ~ 0.4-0.6 keV fearret et alJ 
\2QQ(t . A second Comptonisation component, COMPTT, is also 
needed with ~2.5 keV seed photons coming most likely from 
the NS surface. In our spectra the disk component is too soft 
to be detected and the BMC model interpretes the ~2.5 keV 
seed photons as due to the emission from non-efficiently up- 
scattered seed photons within the bulk flow. The energy sam- 



pling of our data does not allow us to produce a complete 
discussion about the full geometry of the sources and global 
BMC model parameters, beyond the scope of this article, and 
this caveat should be kept in mind in all we discussed above 
that remains, in any case, qualitatively valid. 

4. Discussion 

We have analysed IBIS/ISGRI available public data on a sam- 
ple of twelve persistent LMXBs containing a neutron star. 
We focused our study on the average spectral behaviour of 
the sources and classified them in terms of spectral states. As 
shown by the light-curves we presented, the sources have some 
degree of variability, but the average spectrum is a representa- 
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Fig. 6. Atoll sources: average IBIS/ISGRI spectra and best fit models shown in Tabled 



tion of the sources above 20 ke V and enables the study of the 
global properties. 

4.1. A scenario for spectral evolution in NS LMXBs 

We observe three main spectral states: a very soft state with 
source spectra that can be well described by a single ther- 
mal Comptonisation component, an intermediate state (ther- 
mal Comptonisation plus power-law) and a low/hard state (sin- 
gle power-law). We have successfully studied these three spec- 
tral states in the frame of the Generic Comptonisation Model 
(BMC, TMK96). We present our scenario of the spectral evo- 
lution of NS LMXBs starting from the low/hard state. 



The low/hard state is characterized by a low mass accretion 
rate in the disc. In this case the gravitational energy release 
in the disc is much smaller than the one in the optically thin 
outer boundary of the corona (Compton cloud). The coronal 
outer boundary is p resumably related to the adjustment shock 
dTitarchuk & Fioritoll2004 . The corona completely covers the 
seed photon area (high log A) and the emergent spectrum is a 
result of the up-scattering (Comptonisation) of the seed pho- 
tons in the corona. At the very low level of accretion rate, the 
bulk motion is extremely weak and the emergent spectrum can 
be interpreted as pure thermal Comptonisation with no bulk 
signature (GX 354—0 in Fig. |4}. The bulk motion effect in- 
creases with accretion rate and the power-law signature starts 
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to be visible in the spectrum, the cut-off is moved at a higher 
energy than what expected from pure thermal Comptonisation 
(H 1608-522 in Fig.0}. In either case, one cannot see any trace 
of the seed photons in ISGRI (that give the BB bump). Open 
magnetic field lines are exposed to the observer but the outflow 
is weak because of the low mass accretion rate (there is not 
enough radiation in the disc to launch the wind) and the system 
is radio-quiet. 

In the intermediate state the mass accretion rate increases 
with respect to the low/hard state. It leads to high efficiency of 
the Comptonisation, particularly bulk inflow Comptonisation 
that is seen as an extended hard tail in the spectrum. Thermal 
Comptonisation becomes less efficient because the coronal 
plasma is cooled down by the seed photons coming from the 
disc and NS surface. The corona consists of a quasi-spherical 
component (related to the closed field lines and bulk motion in- 
flow) and of a cylindrical component (related to the open field 
lines and outflow, TS05). The vertical size of the cylindrical 
configuration is suppressed as the mass accretion increases. We 
start to see the seed blackbody bumps (only the higher energy 
one in the IBIS/ISGRI range), because the corona is cooled 
down and becomes more compact. 

The very soft state is characterized by a high mass accre- 
tion rate that is very close to the critical (Eddington) values. 
The emergent spectrum is a sum of two blackbody-like spectra, 
one is related to the Comptonisation of the NS photons (visi- 
ble in the IBIS/ISGRI range) and the other one is related to 
the Comptonisation of the disc photons. The seed photon and 
plasma temperatures differ by a factor of a few. The electron 
plasma and the photons in the Compton cloud are very close 
to equilibrium. In this state the corona is quasi-spherical, there 
is no bulk motion and no radio emission. In fact, the radia- 
tion pressure caused by the strong emission from the NS sur- 
face stops the bulk inflow and the high accretion rate changes 
the configuration of the field lines and the radio emission is 
quenched. 

It is important to note that what plays the key role to sup- 
press the bulk motion is the local radiation pressure, hence the 
local accretion rate impinging the bulk inflow: Sco X-l has a 
total isotropically-estimated luminosity that is higher than e.g. 
GX 3+1, but we speculate that the inflow anisotropy is much 
higher in GX 3+1 than in Sco X-l, resulting in a higher local 
radiation pressure, i.e. bulk suppression. 

4.2. The radio emission - hard X-ray tail connection 

In the scenario proposed in the previous Section there is a 
clear connection between the X-ray spectral states and the radio 
emission: low/hard states are associated to weak radio emis- 
sion, the radio emission increases in the intermediate state and 
then is quenched in the very soft state. 

This trend is clearly met in the observations. We can imag- 
ine a continuous "accretion line" increasing from Atoll sources 
(island state to banana) moving to the Z sources (horizontal, 
normal and flaring branch). Only a few low-dim Atoll sources 
(typical low/hard state) have been de tected in the radio b and 
because of their low radio luminosity JFender & Hendrvl2 000: 



Table 4. Comparison of the average radio (8.5 GHz) and X- 
ray hard tail (40-100keV) flux of the Z and bright Atoll 
sources studied in this p ape r. The radio data are take n from 
iFender & Hendrvl feOOOl) and lMigliari & Fended feOOd) . 



Source 


Radio flux Hard tail flux 
(mJy) 10~ n ergs" 1 cm -2 


Z sources 


Sco X-l 


10 + 3 


18 


GX 17+2 


1.0 + 0.3 


1.9 


GX 349+2 


0.6 + 0.3 


<1.15 


CygX-2 


0.6 + 0.2 


1.7 


GX 5-1 


1.3 + 0.3 


1.9 


GX 340+0 


0.6 + 0.3 


2 


Atolls 


GX9+1 


<0.2 


<0.75 


GX 9+9 


<0.2 


<1.72 


GX 3+1 


<0.3 


<0.43 


GX 13+1 


1.8 + 0.3 


1.74 



1000 



I 00 
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Fig. 7. Average radio luminosity plotted against the average 
hard X-ray tail (40-100keV) luminosity. The correlation is 
clearly visible also if we compare the source fluxes instead of 
the luminosities (see text). 



Migliari & Fender 2006). Z sources in their horizontal branch 
have a spectrum that corresponds to the intermediate state we 
defined here and indee d a clear radio emission has been re- 
ported in this state (Fender & Hendrv 2000; Di Salvo & Stella 
2002, and references therein). The normal and flaring branch of 
Z sources corresponds to the very soft state defined here where 
no radio emission is expected and in fact Z sources in these 
branches have a severely reduced radio emission, if any. 

Z sources are much brighter radio emitters than Atolls and 
G X 13+1 has an important radio emission, thus behaves like a 
Z JFender & Hendrvll2000l iMigliari & Fenderl boQfj). In aver- 
age, unlike bright Atolls, Z sources have a hard tail besides 
the thermal Comptonisation by an optically thick plasma of 
-3 keV JPi Salvo & Stellal2002l and this work); GX 13+1 has 
a hard tail (this work), thus again behaves like a Z. We basi- 
cally see that in GX 13+1, the link radio emission - hard X-ray 
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tail (both "against" its originally declared nature to be an Atoll 
source) reveals itself in a solid way. 

Triggered by the GX 13+1 case, we compared the average 
radio emission of all the sources of our sample with the hard 
X-ray tail that appears on top of the ~3 keV Comptonisation 
spectrum. We note that the hard X-rays from the two soft 
X-ray dim Atoll sources, H 1608-522 and H 1705-440, are 
mainly coming from Comptonisation of soft photons in a hot 
corona (i.e. the overall spectrum is hard, kT e > 40keV, see 
Table |3 and since we cannot disentangle the dominating ther- 
mal Comptonisation component from the dynamical (bulk) 
one, we do not include these two sources in the radio - hard 
tail study. Tab le Elreports the average radio flux as d erived by 
iFender & HendnfeoOOl) and lMigliari & Fended feOOd) and the 
40-100keV flux of the hard tails detected in this work. When 
a hard tail is detected, its 40-100 keV flux is computed af- 
ter fixing the COMPTT normalization to zero. For the sources 
where a hard tail was not detected, we computed 3er upper 
limits (3xerror in the 40-100 keV range), assuming that all 
the flux in the 40-100keV comes from the hard tail. Figure 
is visualising the relation between the average radio and hard 
tail (40-100keV) luminosities computed using the distances 
m Table [I] The correlation between the two is quite strik- 
ing (Spearman rank order significance of 99.6%) and the up- 
per limits are consistent with the trend. We note that we ob- 
tain a very good correlation (Spearman rank order significance 
99.2%) also using the radio and hard tail fluxes instead of the 
luminosities. We chose to plot the luminosity correlation in or- 
der to show the intrinsic properties of the sources. 

In Fig. we have plotted the average radio and X-ray hard 
tail luminosities: the radio luminosities are the superposition of 
optically thick emission (compact jet) and optically thin flar- 
ing activity. The X-ray hard tail luminosities are the average of 
all hard tail states (from maximum to minimum strength). We 
are aware that these correlations are based on data that are not 
taken simultaneously, but they are all based on average fluxes 
that we can consider a good representation of the physics in- 
volved in these objects. For the sources where a hard tail has 
been observed, we compare the radio emission with what is 
left in the X-ray spectra of the sources once we remove the 
dominant Comptonisation component. In this respect, our de- 
tection of a hard tail in the spectrum of GX 13+1 is important: 
GX 13+1 has a similar X-ray flux (and luminosity) as the re- 
maining bright Atoll sources (GX 3+1, GX 9+1, GX 9+9, see 
TableQ but a much higher average radio and hard tail emission 
(Table[4}, re-confirming the radio hard X-ray tail connection. 

In the presence of a cut-off, the hard tail can be explained 
via Comptonisation of soft seed photons in the jet and/or in the 
corona. For the case where no cut-off is detected, many models 
have been propo sed: Com ptonisation by a hybrid (thermal non- 
thermal) corona (Coppi 1999), synchrotron emission from the 
electrons of the jet ( Markoff et al. 2005) or bulk motion inflow 
(dynamical) Comptonisation (TMK96, used in this work). 

The real test of any of these models can be done using the 
(variability) analysis of the power density spectrum (PDS) of 
the hard tail emission. The characteristic (break and QPO) fre- 
quencies of PDS do determine the geometric size of the con- 
figuration where the hard tail emission is formed. TS05 made 



this type of analysis for 4U 1728-34 and found evidence that 
the hard tails are formed in the compact Compton cloud with 
geometry changing from cylindrical-like in the low/hard state 
to the quasi-spherical one in the high/soft state. Note that the 
X-ray-radio correlation along with the QPO-radio correlation 
is well established in 4U 1728-34. TS05 presented an expla- 
nation of these correlations in the framework of an oscillation 
model using the observed correlations of QPO low frequencies 
and their ratio. 

The correlation we find between the radio and hard X-ray 
tail emission suggests that the hard tail formation area and 
the source of energetic electrons, ultimately causing the radio 
emission, are closely connected. The most probable site of this 
configuration is the NS magnetosphere. It can be suggested that 
the open magnetic field lines of the NS magnetosphere are the 
base of the jet seen in the radio emission. An increasing ac- 
cretion rate leads to a more efficient radio emission (low/hard 
state to intermediate state) up to a point where the extremely 
high accretion rate (very soft state) changes the configuration 
of the field lines and the radio emission is quenched. 
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